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Abstract. We investigate the ionization and the opacity of warm, dense helium under conditions found in 
the atmospheres of cool white dwarf stars. Our particular interest is in densities up to 3 g/cm 3 and tempera- 
tures from 1000K to 10000K. For these physical conditions various approaches for modeling the ionization 
equilibrium predict ionization fractions that differ by orders of magnitudes. Furthermore, estimates of 
the density at which helium pressure-ionizes vary from 0.3 to 14 g/cm 3 . In this context, the value of the 
electron-atom inverse bremsstrahlung absorption is highly uncertain. We present new results obtained from 
a non-ideal chemical model for the ionization equilibrium, from Quantum Molecular Dynamics (QMD) 
simulations, and from the analysis of experimental data to better understand the ionization fraction in fluid 
helium in the weak ionization limit. 
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INTRODUCTION 

We are interested in the opacity of He and He/H mix- 
tures for densities up to 3 g/cm 3 and temperatures of 
T = 1000K to 10000K to model the atmospheres of 
very cool white dwarf stars. Our work is motivated 
by the fact that there exists a wide range of predic- 
tions for the density 

TABLE 1. Density ppj for the pressure ion- 
ization of helium from experiments (exp), 
chemical models (chem) and a T = K quan- 
tum mechanical calculation (LMTO) 



Refference ppj (g/cm 3 ) Source 
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chem 
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chem 
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chem 
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>6.5 


chem 
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7.5-14 


LMTO 


[1] 


~ 1.5-2 


chem & exp 



at which helium pressure-ionizes (Table 1). This 
translates into large uncertainties in the number of 
free electrons and the value of the free-free opacity in 
the regime of our interest. We approach this problem 
by constructing more reliable models for ionization 
equilibrium of He on the basis of a chemical model, 
QMD simulations, and available experimental data 
[1]. 



THE CURRENT STATE OF MODELING 

There is a large span in predictions for the density at 
which the pressure ionization phenomenon occurs in 
helium (Table 1.). 

For the chemical models [l]-[5], this arises from 
oversimplified treatments of the physics and the pres- 
ence of free parameters. The interactions between 
atoms and ions are usually described by a polariza- 
tion potential with a hard sphere cut off at short dis- 



tance (e.g. [5]), which varies between models. Often 
they are simply neglected (see [4]). In particular, It 
is not correct to treat the e — He interaction with a 
polarization potential [7], as we discuss below. 

The more sophisticated quantum calculations of 
the r = 0K solid predict a high density for pressure 
ionization [6], which we confirm with QMD simula- 
tions. On the other hand the only published exper- 
imental data on helium for the regime of pressure 
ionization [1] suggest that this phenomenon occurs 
at density ~ 1.5 -2.0g/cm 3 . 



THE ELECTRON-HELIUM INTERACTION 

The energy of a free electron in cryogenic helium 
has been shown experimentally to be as large as 
V e <~ +1.5 eV (positive) at densities of ~ 0.2 g/cm 3 
[8]. This is in the quantitative agreement with the 
prediction of the Lenz potential energy [9] 

V e = 2%h 2 n{fte)a/m e , (1) 

where a ~ 1.3 a.u. is the scattering length [9], and 
n(He) is the number density of helium. This is be- 
cause the e He interaction is strongly repulsive for 
r > 1 a.u. [10], and significantly different from the 
polarization interaction [7] (Fig. la, solid & dotted 
curves). The quantum mechanical calculations for 
the energy of an electron in dense helium, using the 
potential of [10], are in good agreement with exper- 
imental data [11]. This shows that free electrons in- 
teract strongly with the atoms in dense He, but not 
through an attractive polarization potential. This is 
an important element in chemical models. 

To introduce the interaction of quantum mechani- 
cal electrons with the atoms in a chemical model of 
the thermodynamics, we performed Density Func- 
tional Theory (DFT) calculation of the free energy 
of free electron in dense He. The results are well ap- 
proximated by the Lenz formula (Fig. lb), which we 
adopt in the chemical model. 

THE CHEMICAL MODEL FOR THE 
IONIZATION EQUILIBRIUM 

We constructed a chemical model considering the 
following species: He,He + ,HeJ,e . The interac- 
tions between atoms and ions are described by the 



potentials shown in Fig. la. In the chemical picture, 
the non-ideal effects can be treated as a shift in the 
ionization/dissociation energies [3]. Considering the 
reactions He <-> He + + e~ and He \ <-► He + He + the 
shifts are Al x = J u m ' </ (He) - n nid (He+) - \i nid {e) and 
AI 2 = M mW (HeJ) - jU m ' d (He) - jU m ' d (He+). The non- 
ideal contributions to the chemical potentials for the 
He atom and the ions where obtained through the 
numerical solution of the Ornstein-Zernike equation 
in the Percus-Yevick approximation [15]. This ap- 
proach is valid in the low ionization limit. The re- 
sults are shown in Fig. 2a, where they are compared 
with the shift in the ionization energy extracted from 
conductivity data [1] through the relation O exp / a id = 
gMi/kaT > w here is the ideal gas conductivity. The 
agreement with the experiment is good. In this model 
He pressure-ionizes at ~ 2 g/cm 3 because of strong 
attractive He — He + interaction (Fig. la), which fa- 
vors ionization. This model is questionable, however, 
because this He — He + potential inevitably leads to a 
bound state (HeJ ), which would greatly reduce the 
interaction with the other neighboring He atoms. 



QUANTUM MOLECULAR DYNAMICS 
RESULTS 

We also conducted QMD-DFT calculations of dense 
helium, using the Viena Ab-initio simulation pack- 
age. The results in terms of band gap (ionization 
energy) is presented on Figure 2b. We find that in 
this model, helium pressure-ionizes above a den- 
sity of 10 g/cm 3 . The resulting conductivities at the 
experimental conditions are 2 orders of magnitude 
smaller and show a very strong temperature depen- 
dence rather than the strong density dependence re- 
ported in [1]. It is well known that band gaps are un- 
derestimated using GGA functional. More accurate 
functionals would only increase the gap and worsen 
the disagreement with the experimental data. 



THE FREE-FREE OPACITY FROM HELIUM 

The free-free (inverse bremsstrahlung) opacity of 
dense helium is determined by the ionization frac- 
tion and e — He collisions. In a dense medium, the 
collisions can be described with the classical Drude 
model [17]. On Figure 3 we compare the QMD ab- 
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FIGURE 1. (a) The interaction potentials for the chemical model. The lines represent the following potential curves: e — He 
of [10] (solid) and [7] (dotted), He -He of [12] (dashed), He-He+ of [13] (long dashed), andHe-He+ of [14] (dash-dotted), 
(b) The DFT free energy of a free electron in helium. The line represents the Lenz potential energy (Eq. 1). 




FIGURE 2. (a) The change in the ionization/dissociation potential predicted by the chemical model. The lines represent the 
results for A/j (solid, dotted) and A/2 (dashed, dash-dotted) for T = 1.0 eV and T = 2.0 eV respectively. The filled circles 
represent the value of A/j A/ ' extracted from the experimental data of [1], (b) The QMD band gap for T = 0.5 eV (solid line) 
and T = 1.5 eV (dotted line). The filled circles represent the effective value of ionization energy 24.6 eV— AI* xp . 
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FIGURE 3. The absorption coefficient of He obtained 
from QMD simulations (solid lines) for temperature T = 
0.5 eV and densities p = 0.5,2.0 and 4.0 g/cm 3 (from the 
bottom to top). The dotted line represents the standard free- 
free absorption coefficient for p = 0.5 g/cm 3 [19], adopt- 
ing the density of free electrons of the QMD simulation. 

sorption coefficient [18] with the standard free-free 
absorption coefficient corrected for the correlations 
in dense fluid but not the e — He collision frequency 
[19]. For the astrophysical applications, we are in- 
terested in the spectral region with photon energies 
from zero to 4 eV. The standard free-free frequency 
behaviour is erroreous, because the absorption pro- 
cesses for photons with small energies are driven by 
the frequent e — He collisions rather than the slowly 
varying electric field of radiation. 



CONCLUSIONS 

Our goal is to calculate the free-free absorption from 
dense, non-ideal, weakly ionized helium. This source 
of opacity depends on the ionization fraction and the 
frequency of e He collisions. We use two differ- 
ent approaches to solve for the number density of 
free electrons in dense helium: a chemical model and 
QMD-DFT simulations. In the first method we have 
introduced a new description of the e — He interac- 
tion and use ab initio potentials for the He-ions in- 
teractions. This gives us a good agreement with the 
experiment, but the result is driven by a questionable 
He — He + potential. On the other hand, the more so- 
phisticated QMD simulations suggest that pressure 
ionization of He does not occur below 10 g/cm 3 , 
which is inconsistent with the conductivity data. This 
large discrepancy between models and the conduc- 



tivity measurements is astrophysically significant. 
New experiments that probe the pressure ionization 
of dense He are essential in finding a resolution of 
this problem. 
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